Statistics of active vs. passive advections in magnetohydrodynamic turbulence 
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Active turbulent advection is considered in the context of magneto- hydrodynamics. In this case, 
an auxiliary passive field bears no apparent connection to the active field. The scaling properties 
of the two fields are different. In the framework of a shell model, we show that the two-point 
structure function of the passive field has a unique zero mode, characterizing the scaling of this field 
only. In other words, the existence of statistical invariants for the decaying passive field carries no 
information on the scaling properties of the active field. 
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In the context of turbulent advection, the understand- 
ing of fluid turbulence has greatly improved in the re- 
cent years QJ . The anomalous scaling has been shown to 
be universal and connected to the existence of statisti- 
cal integrals of motion 0] . In , it was shown that the 
statistically conserved structures of decaying passive tur- 
bulence dominate the statistics of forced turbulence, thus 
offering a rather general framework for understanding the 
universality of anomalous scaling in forced turbulence. 

Let <j) be a decaying field transported by a stationary 
turbulent flow. The linearity of the advection implies the 
following relation for the correlation functions : 

^{r U t)...(t>{r N ,t)) = 

dqPW(r,t\q,to)(<t>(qi,to)...<fi(qN,t )) , (1) 



where we used the compact notation r = r\ , . . . , rjv to 
denote a collection of N position vectors. Equation 
tells us there exists a linear operator J>( N ) that propa- 
gates the nth order correlation function from time to to 
time t. Without fresh input, that is in the absence of 
forcing, the correlation functions of <f> decay due to dissi- 
pative effects. Nevertheless, as conjectured in [3(, there 
exist special functions Z^ N > that are left eigenfunctions 
of eigenvalue 1 of the operator V^ N \ 



ZW(r)= dqPW(q,t\£,t )ZW(q), 



such that 



(2) 



drZ^ N \f){cj>{q l ,t)... ( j>{q N ,t)) (3) 



is preserved in time. I^ N ' and Z^ N > are respectively 
called a statistical integral of motion and a statistically 



preserved structure of order N, also referred to as zero 
modes [TEj . 

Now, consider the same passive advection problem 
with an external forcing, such that the system reaches 
a stationary state. Define, the correlation function of 4> 
in that stationary state to be, 



(4) 



thomas.gilbert@inln.cnrs.fr 



' dhruba@physics.iisc.ernet.in 



where the symbol (■)/ denotes averaging over the sta- 
tistical stationary state. It was conjectured in Ref 
that the anomalous part of F^ N \f) is dominated by 
the leading zero modes of the decaying problem, i. e. 
Z W ~FW. The conjecture was verified in the context 
of a shell model for passive scalar advection. 

In subsequent studies, it was discovered that the ex- 
istence of statistical invariants of the motion for passive 
turbulence may help understand the statistics of active 
turbulence, a case where the advected quantity affects 
the dynamics of the advecting field. In Q and 0, El, 
the case of thermal convection in the Boussinesq approx- 
imation was studied. There it was shown that the scal- 
ing of the active field is also dominated by the statisti- 
cally preserved structures of auxiliary passive fields. It 
is yet unclear how general this connection between the 
statistics of active and auxiliary passive fields is. The 
case of 2-dimensional magneto-hydrodynamics is reveal- 
ing. Indeed, in this case active and passive fields have 
very different scaling behaviors. As shown in 0, the 
transported fields cascade in different directions. It is 
therefore unexpected that the statistics of the auxiliary 
passive field holds information on the statistics of the ac- 
tive field. The two fields have different scaling properties. 
In [fj, despite this difference, the claim was made that 
the analogy does hold in the sense that there exist sub- 
leading zero modes of the propagator of the correlation 
functions of the auxiliary passive field with the scaling of 
the correlation functions of the active field. 

The purpose of this note is to show that this is actually 
not the case. To this end, we will limit our investigation 
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to the case of the second order structure functions of a 
shell model of 2-dimensional magneto-hydrodynamic tur- 
bulence. The case of the second order structure function 
is the simplest one. Because of the absence of geome- 
try, its scaling is non-anomalous 0, 0, , which in the 
language of zero modes implies non-degeneracy, that is 
there is a unique conserved structure associated to the 
two-point statistical invariant. The same holds in the 
language of shell models, where the only two point func- 
tion of a passive scalar field 6 n associated to the scale k n 
is (\9 n \ 2 ), in contrast to higher order structure functions, 
e. g. the fourth order for which we have contributions 
from (|#„| 4 ), (|#„-i| 2 |#n+i| 2 ), etc. Using methods simi- 
lar to those used in pjj, we will construct the operator 
propagating the second order structure functions, and 
will demonstrate that the auxiliary passive field has the 
same statistical integral of motion as other shell mod- 
els of passive advection. It will be inferred that there is 
no sub-leading zero-mode with the scaling of the (active) 
magnetic field. 

In analogy to other models 0, 0, the follow- 
ing two sets of equations generalize the usual Sabra shell 
model 14] for the turbulent velocity field to magneto- 
hydrodynamic turbulence (we omit dissipative terms) : 



du n 
dt 



db n 
dt 



i[k n +l{K+l u n+2 ~ K+l h n+2) 

-(e + l)fc„«_ 1 u„ + i - b* n _ x b n+x ) (5) 
-ek n -i(u n -2U n -i - K-2K-1)] + fn, 

i[-(e + S)k n+1 (u* n+1 b n+2 - b* n+1 u n+2 ) 
+5k n (u* n _ l b n+ i - bn^Un+i) (6) 
+ (<5 - l)fc„_i(u„_ 2 &n-l - 6„-2«n-l)] + fn- 

Here u n (the velocity field) and b n (the magnetic field) 
are complex variables defined on a discrete set of shells 
indexed by the integer n whose associated wave-number 
k n = ko\ n , A > 1 (hereafter taken to be A = 2). /„ and 
f' n are two forcing terms which are taken to be stochastic 
white noises with identical statistics and concentrated on 
a limited number of neighboring shells (n = 5, 6, 7 in our 
numerical experiments). The model's parameters e and S 
are conveniently parametrized in the following way. The 
three-dimensional model for which e < reads (a > 0) 



e=-\- a , 
5= (l + A")- 1 



(3D model) 



(7) 



The two-dimensional problem on the other hand has e > 
and reads 



e = A" Q 
S = -(\ c 



1) 



(2D model) 



(8) 



Correspondingly, we have the following quadratic dynam- 
ical invariants (i. e. time-invariant in the limit of zero 
viscosity and zero external forcing). 



E 



(total energy), 



(9) 



K = ^K«&„) (cross helicity), (10) 

n 

H = ^sign(<5)"fc- a |&„| 2 (magnetic helicity) .(11) 



Thus, in both two- and three-dimensional models, the 
equations JSJ EJ) have one single free parameter, a > 0. 
The two- and three-dimensional models actually have 
very different dynamical behaviors, see ^]J. Only the 
two-dimensional model can sustain a stationary state and 
we will limit ourselves to this case. As dimensional anal- 
ysis shows [g , the conservation of the first two invariants 
implies that u n and b n must both have Kolmogorov scal- 
ings, (|u„| 2 ), (\b n \ 2 ) ~ fc,7 2 ^ 3 , for which the corresponding 
fluxes are constant. This is indeed what has been mea- 
sured for similar models [ll|, where both fields appear 
to display the same anomalies. Further, as shown in |(|, 
the conservation of the third invariant allows for another 



-2/3 



for which the magnetic helicity 



scaling, (|6„| 2 ) 
flux is constant. But since this scaling is incompatible 
with the conservation of the two other invariants, it is 
not relevant to the statistics of the magnetic field. 

However if one considers a passive auxiliary field obey- 
ing an equation identical to Eq. JBJ for the 2-dimensional 
case, 



da n 
~dT 



i[X "(A"-l) 1 k n+1 (u* n+1 a n+2 - a* l+1 u n+2 ) 
~ (A Q - l) _1 fc„«_ 1 a„ + i - a^tin+i) 

-A Q (A Q - iy 1 k n - 1 (u n - 2 an-i - a n _ 2 tt7i-i)] , 

(12) 



the only relevant invariant is the equivalent of the mag- 
netic helicity Eq. l|llfl . to which a constant flux is asso- 
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is 



ciated. Thus the dimensional scaling (|a„| 2 ) ~ kr, 
expected to be observed. And that is indeed what was 
found in 0. 

Notice though that the linearity of Eq. i|12|) allows for 
the substitution VVt = A _ "™/ 2 a„. With this new variable, 
Eq. l(T2"|) takes the form, 



dt 



n — 2 u n—l ) 

+ B(k n+1 ll)* n+1 U n+2 - fenVC-l^n+l) 

+ C{k n U* n _ l 1p n+ l + k„-lU n ^2lpn-l)}, (13) 



which describes the advection of a scalar for which the 
quadratic invariant is |"!/vi| 2 , which is similar to the 
shell model with only nearest neighbor interaction con- 
sidered in 0, . The coefficients in Eq. l(T2")l correspond 
to the choice A = (A" — 1)~\ B = -A~ Q / 2 (A a - l)-\ 
andC= -A Q / 2 (A Q -1)" 1 . This shows that as far as shell 
models are concerned, the difference between a passively 
advected vector and passively advected scalar is just a 
numerical factor. Henceforth we do not make a distinc- 
tion between the two and use the name "passive field" 
for both. 
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Consider now the equivalent of Eq. for the propa- 
gation of the second order structure functions (|# n | 2 ) in 
the decaying problem - including dissipative terms on the 
RHS of Eq. ((13(1 . Following notations similar to those 
used in [lOj , we can write the equation of motion for ip n 
under the form 



dt 



£"n,m^m : (1^0 

with the solution 

ipn(t) = T+ ^ exp / dsC(s) I tp m (to), 

v L^io J J n,m 

= TZ„^ m {t\to)lp m (to), (15) 

(T + denotes the time ordering operator). Letting 

v^ m (t\t ) = (Hn, m (t\to)K, m (t\to)) > (is) 

the propagation of second order structure functions obeys 
the following equation : 

(\Mt)\ 2 ) = Y, V n%(t\to)W m (ta)\ 2 ) ■ (17) 

m 

The form of the operator was discussed in It is 
a matrix whose elements can be obtained by propagating 
an initial condition concentrated at a given shell. Simi- 
lar considerations hold for the models considered here. In 
Fig. \I\we show these elements for successive times, start- 
ing from an initial conditions at shell 20. The model we 
used is Eq. (|13|) for the parameters corresponding to Eq. 
(|12fl and advected by the magneto-hydrodynamic fields 
Eqs. (JSJEl). The parameters of the simulation are given 
in the figure caption. 

In analogy to Q, the statistical invariant for the 
passive scalar field is 



i^(t) = J2z^{\Mt)\ 2 ) 



(18) 



where Z n ' is a left-eigenfunctions of the operator V^ 2 \ 
with the scaling of the second order structure function of 
the forced problem. The invariance of I^ 2 ' is most eas- 
ily demonstrated by re-scaling the decaying second order 
objects according to Eq. 1)18(1 . The curves indeed col- 
lapse if the ordinate is shifted with the appropriate time 
dependence. This is shown in Fig. |3 and is analogous 
to Fig. 4 in [ijj. Thus in the language of the passive 
magneto-hydrodynamic model Eq. 1(12(1 . we have that 



£<KlV fc n(Kr(*)>A" 



(19) 



is a statistical integral of motion. To claim that the scal- 
ing of the magnetic field structure function (\b n \ 2 ) is a 
sub-leading zero mode of is equivalent to claiming 
that 



J2(\bn\ 2 )/K(\a n \ 2 (t))/k« 



(20) 




FIG. 1: The elements of V^ 20 for the scalar field Eq. 
advected by the MHD field Eqs. (0 [§J , where the parameter 
was chosen to be a — 2. The times displayed are .4, .55, .7, 
.85 and 1 respectively (measured in the natural time units 
of the model). The horizontal axis corresponds to the shell 
numbers. The units on the vertical scale are arbitrary. The 
simulation was done using a total of 35 shells, with the first 
shell wave-number ko = 1/16 and a shell spacing of A = 2. 
All the fields were dissipated on the small scales with a term 
vkn, with v — 10 -12 . The advecting fields were forced on 
shells 5 — 7 with white delta correlated noise of amplitudes 
I/a/2, l/2\/2 and 1/4 respectively. Moreover the phase of 
the forcing on shell 7 was taken to be equal to the sum of the 
phases of the forcings on shells 5 and 6. 
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FIG. 2: The curves of Fig. Qcollapsed according to Eq. 1(18^ . 



is a statistical integral of motion. Clearly, in view of 
Fig. El this cannot be the case and the collapse will not 
occur should the scaling exponent of be replaced by 
another one In this line of thought, it is perhaps 

worthwhile pointing out that the collapse as seen in Fig. 
El would not be possible should there be zero- modes with 
distinct scaling exponents. Indeed, as seen from Eq. l(T%)l . 
the collapse occurs provided {\ip n (t)\ 2 ) "falls" precisely on 



the right eigenmode with a scaling identical to Zn 



(2) 
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To conclude we emphasize that the linearity of the pas- 
sive advection models assigns them to a narrow class of 
equivalence. The passive magneto-hydrodynamic model 
is in fact equivalent to a scalar advection model for which 
the statistical invariants have already been investigated 
in some detail [TfT ]. In magneto- hydrodynamics, active 
and passive fields have different scaling properties. The 
arguments that were used in the framework of thermal 
convection to account for the anomalous scaling of the ac- 
tive field in terms of a passive auxiliary one do not carry 
over to magneto-hydrodynamics. The claim that one can 
nevertheless account for the scaling of the (active) mag- 
netic field by considering sub-leading zero-modes of the 
operators propagating the decaying correlation functions 
was proven wrong. In view of the form of the second 
order propagator, it is clear that there are no zero mode 
but the one whose scaling is that of the passive field. 
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